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Abstract. We show that (n, 2") additive codes over GF(4) can be repre- 
sented as directed graphs. This generalizes earlier results on self-dual ad- 
ditive codes over GF(4), which correspond to undirected graphs. Graph 
representation reduces the complexity of code classification, and enables 
us to classify additive (n, 2") codes over GF(4) of length up to 7. From 
this we also derive classifications of isodual and formally self-dual codes. 
We introduce new constructions of circulant and bordered circulant di- 
rected graph codes, and show that these codes will always be isodual. A 
computer search of all such codes of length up to 26 reveals that these 
constructions produce many codes of high minimum distance. In par- 
ticular, we find new near-extremal formally self-dual codes of length 11 
and 13, and isodual codes of length 24, 25, and 26 with better minimum 
distance than the best known self-dual codes. 

Keywords: additive codes; quaternary codes; classification; graphs; cir- 
culant codes; formally self-dual codes 



1 Introduction 

An additive code, C, over GF(4) of length n is an additive subgroup of GF(4) n . 
We denote GF(4) = {0, 1, u, to 2 }, where J 1 = u + 1. C contains 2 k codewords for 
some < k < 2n, and can be defined by a k x n generator matrix, with entries 
from GF(4), whose rows span C additivcly. C is called an (n,2 k ) code. In this 
paper we will only consider (n,2 n ), or half-rate, codes. 

The Hamming weight of u G GF(4) ra , denoted wt(tt), is the number of 
nonzero components of u. The Hamming distance between u and v is wt(u — v). 
The minimum distance of the code C is the minimal Hamming distance between 
any two distinct codewords of C. Since C is an additive code, the minimum dis- 
tance is also given by the smallest nonzero weight of any codeword in C. A code 
with minimum distance d is called an (n, 2 k , d) code. The weight distribution of 
the code C is the sequence (A ,Ai, . . . ,A n ), where Ai is the number of codewords 
of weight i. The weight enumerator of C is the polynomial 

n 

w c {x,y)^Y. A ^ n ~ l y l 

i=0 
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Two additive codes over GF(4), C and C , are equivalent [T] if and only if the 
codewords of C can be mapped onto the codewords of C by a map that consists 
of a permutation of coordinates (or columns of the generator matrix), followed 
by multiplication of coordinates by nonzero elements from GF(4), followed by 
possible conjugation of coordinates. Conjugation of x G GF(4) is defined by 
x = x 2 . For a code of length n, there is a total of 6 n n! such maps. The 6 
possible transformations given by scaling and conjugation of a coordinate are 
equivalent to the 6 permutations of the elements {1,uj,ui 2 } in the coordinate. 
A map that maps a code to itself is called an automorphism of the code. All 
automorphisms of C make up the automorphism group, denoted Aut(C). We can 
use the computational algebra system Magma 2 to find the automorphism group 
of a code. Since, at this time, Magma has no explicit function for calculating the 
automorphism group of an additive code, we use the following method, described 
by Calderbank et al. [3]. We map the (n, 2 k ) additive code C over GF(4) to the 
[3ra, k] binary linear code /3(C) by applying the map i— » 000, 1 i— > Oil, u> \— > 101, 
uj 2 >— > 110 to each generator of C. We then use Magma to find Aut(/3(C)) n 
Aut(/3(GF(4)™)), which will be isomorphic to Aut(C). 

The trace map, Tr : GF(4) -» GF(2), is defined by Tr(x) = x + x. The 
Hermitian trace inner product of two vectors over GF(4) of length n, u = 
(ui,u 2 , . . . ,u n ) and v = (vi,v 2 , . . - ,v n ), is given by 

n n 

u*v = Tr(tt • 17) = Tr(uiVi) — '^^(u i vf + ufvi) (mod 2). 
»=i i=i 

We define the dual of the code C with respect to the Hermitian trace inner 
product, C 1 - = {u e GF(4) n | u * c = for all c G C}. C is self-dual if C = C ± , 
formally self-dual [4] if We(x,y) = W c ±{x,y), and isodual if C is equivalent to 
C . All self-dual codes are isodual, all isodual codes are formally self-dual, and all 
formally self-dual codes are half-rate codes. The set of linear half-rate codes over 
GF(4) is a small subset of the additive half-rate codes of even length. Optimal 
linear half-rate codes over GF(4) of length up to 18 were classified by Gulliver, 
Ostergard, and Senkevitch [5]. The set of half-rate additive codes contains all 
self-dual, isodual, and formally self-dual additive codes, as well as all half-rate 
linear codes. 

It follows from the Singleton bound [6] that any half-rate additive code over 
GF(4) must satisfy 



d < 




C is called extremal if it attains the minimum distance d given by the Singleton 
bound, and near-extremal if it has minimum distance d—1. If a code has highest 
possible minimum distance, but is not extremal, it is called optimal. Han and 
Kim |4l7j showed that there are no extremal formally self-dual codes of length 
n > 8, and no near-extremal formally self-dual codes of length n = 16, n — 18, 
or n > 20. Tighter bounds on the minimum distance of self-dual additive codes 
over GF(4) were given by Rains and Sloane [HI Theorem 33]. 
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One of the motivations for studying self-dual additive codes over GF(4) has 
been the connection to quantum error- correcting codes [3j. Non-self-dual additive 
codes cannot be applied as quantum codes in the same way, but are interesting 
for other reasons. Han and Kim [4l7j studied formally self-dual additive codes 
over GF(4), and showed that some of these codes have higher minimum distance 
than the best self-dual codes of the same length. Additive codes may also be 
better than the best linear codes of the same length. It is known that some 
strong binary codes can be projected onto additive codes over GF(4) A 
connection between formally self-dual codes over GF(4) and lattices has also 
been shown [TIT]. We have previously studied the connection between self-dual 
additive codes over GF(q 2 ), for any prime power q, and weighted graphs 
Such generalizations could also be considered for additive codes in general, but 
will not be discussed in this paper. 

Let t n be the number of incquivalent codes of length n. To find one code 
from each of the t n equivalence classes, i.e., to classify the codes of length n, is 
a hard problem. We have previously classified all self-dual additive codes over 
GF(4) of length up to 12 [TJ], by using the fact that all such codes can be 
represented as undirected graphs |13ll4ll5ll6j , and that an operation called local 
complementation (LC) generates orbits of graphs that correspond to equivalence 
classes of codes |13ll6j . 

The main result of this paper is to show that additive (n, 2 n ) codes over 
GF(4), except for some special cases, have representations as directed graphs. 
This basically transforms the problem of classifying such quaternary codes to 
a binary problem, with reduced complexity. We show that an algorithm by 
Ostergard [17] for checking equivalence of linear codes also works for additive 
codes over GF(4). By using this algorithm, and the fact that codes correspond 
to directed graphs, we are able to classify additive (n, 2 n ) codes over GF(4) of 
length up to 7. We find that there are more than two million inequivalent codes 
of length 7. We have also checked which codes are formally self-dual, isodual, 
or self-dual, and give details of this enumeration. We introduce circulant and 
bordered circulant directed graph codes, and a computer search of all such codes 
up to length 26 reveals this subclass of additive half-rate codes to contain many 
codes with high minimum distance. Due to the structure of the generator ma- 
trices, codes from these constructions will always be isodual, and hence also 
formally self-dual. We construct new near-extremal formally self-dual codes of 
length 11 and 13, which were previously unknown [4]. This also answers in the 
affirmative the open question of the existence of an additive (13, 2 13 , 6) code [18]. 
Finally, we find isodual codes of length 24, 25, and 26 with minimum distance 9. 
The best known self-dual codes of these lengths have minimum distance 8. 

2 Directed Graph Representation 

A directed graph is a pair G = (V, E) where V is a set of vertices, and E C V x V 
is a set of ordered pairs called edges. A graph with n vertices can be represented 
by annxn adjacency matrix r, where 7y = 1 if (i,j) <E E, i.e., if there is a 
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directed edge from i to j, and jij = otherwise. We will only consider simple 
graphs, where all diagonal elements of the adjacency matrix are 0. The special 
case where we always have an edge (j, i) G E whenever there is an edge (i, j) G E, 
i.e., the adjacency matrix is symmetric, is called an undirected graph. The in- 
neighbourhood of v G V, denoted NI V C V, is the set of vertices i such that 
there is a directed edge (i, v) G E. Similarly, NO v C V is the out-neighbourhood 
of v, i.e., the set of vertices i such that there is a directed edge (v, i) in E. \NI V 
is the indegree of v, and \NO v \ is the outdegree of v. Two graphs G = (V, E) 
and G' = (V, E') arc isomorphic if and only if there exists a permutation tt of 
V such that (w, v) G -E (7t(m), 7r(w)) G -E'. A directed graph is connected, 

(also known as weakly connected), if we can reach any vertex starting from any 
other vertex by traversing edges in some direction, i.e., not necessarily in the 
direction they point. 

Definition 1. A directed graph code is an additive (n, 2") code over GF(4) that 
has a generator matrix of the form C = r + loI, where r is the adjacency matrix 
of a simple directed graph and I is the identity matrix. 

Proposition 1. Given a directed graph code C with generator matrix C = T + 
cul , its dual C 1 - is generated by C T . 

Proof. We must show that for any c G C and any d G C ± , the trace inner 
product c * d = Tr(c • d) = 0. Let c = aC and d = bC T , with a, b G GF(2) n . 

Then c*c' = (aC)*(bC T ) = Tr((aC)_(feC^]) = Tr {{aC){bC~ T ) T ) = Tr{aCCb T ), 
which must be if all elements of CC are from GF(2). This is clearly the case, 
since CC = (r + Lul)(r + lu 2 I) = T 2 + T + I. □ 

Theorem 1. Given an additive (n, 2") code C over GF(4) with generator matrix 
C , there exists a directed graph code equivalent to C, except when C = A + coB is 
such that all 2™ sets of n columns ({a\, bi}, {a 2 , b 2 }, ■ ■ ■ , {a„, b n }) of (A \ B), 
where {a^bj} means that we choose either di or bi, are linearly dependent. 

Proof. We can write C = A + ojB, with (A | B) a binary nx2n matrix. From the 
fact that the rows of C additively span a vector space of dimension n, it follows 
that (A | B) has full rank. If the n x n submatrix B also has full rank over 
GF(2), we simply perform the basis change B~ X (A | B) = (r'\T). Any non-zero 
elements on the diagonal of i~" can simply be set to zero, effected by conjugating 
the corresponding coordinates of J" +coI, to obtain an equivalent directed graph 
code generated by r + u)I. 

In the case that B does not have full rank, we must show that there is a code 
C, equivalent to C, generated by A' + uuB' where B' does have full rank. Then 
we can apply the method described in the first part of this proof to obtain the 
graph form. Let the columns of A be denoted (oi, ■ ■ ■ , a n ) and the columns 
of B be denoted (6i, 62, • • • , b n ). Observe that multiplying column i of C by w 2 , 
followed by conjugation of the same column, has the effect of swapping columns 
di and bi in ( A \ B) . Since we know that at least one of the 2" possible choices 
({ai, t>i}, {a,2, 62}, ■ • • i {a-m b n }) span a vector space of dimension n, we can find 
a matrix (A' \ B') where B 1 has full rank. □ 
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It follows from Prop. [T] that a directed graph code is self-dual if and only if its 
generator matrix is symmetric, i.e., it is in fact an undirected graph code. The 
fact that all self-dual additive codes over GF(4) can be represented as undirected 
graphs is well known [13114115116] , and was used to classify all self-dual additive 
codes up to length 12 [TJ]. Theorem [1] is a generalization of this result to the 
much larger classes of directed graphs and half-rate additive codes over GF(4). 
As stated in Theorem [TJ some special codes do not have graph representations. 
These are codes that will typically not be of interest, such as codes with a 
generator matrix that contains an all-zero column or a set of linearly dependent 
binary columns (up to scaling by uj or uj 2 ). 

Example 1. We consider an additive (7, 2 7 ,4) code, C, generated by 



C 



(uj 2 

1 



1 
1 



(A\B) = 



10 

1 

uj uj 2 

1 1 uj 2 

UJ 1 LU 2 

1 1 1 UJ 2 

1 uj 1 



uj\ 
1 



UJ UJ 

uj uj 2 
uj 2 



A + ujB, 



/l 1 
10 111 
110 
1110 1 
10 1110 

1111111 

\0 1 1 1 1 



1/ 

1 1 l\ 

10 111 
111 
10 110 
10 
1 1 0/ 



We swap column ae with be and column aj with bj to get the matrix 



(A' | B') 



^10 


1 





1 1 


1 0\ 


1 


10 


11 








1 


1 1 


10 110 


1 


1 


1 


1 1 


10 1 


1 


1 


1 


1 


10 110 


1 1 1 


1 


1 





111 


[p 1 





1 


1 


1 1 1/ 



(J" | 7) = B'^iA' I B') = 



obtain the matrix 








A o o l o l l 


10 










10 10 1 


10 











10 1110 


10 











10 10 10 


1 











1110 





1 








110 11 








1 





\0 1 1 1 1 











V 



() 



By setting the diagonal of r' to zero, we get the adjacency matrix of a simple 
directed graph, 

/0 1 1 l\ 
10 10 1 
10 110 

r = ioiooio 

1110 
110 1 
\0 1 1 1 0/ 

This graph is shown in Fig. [1] r + u>I generates a (7, 2 7 , 4) directed graph code 
equivalent to C. 



Fig. 1: Directed Graph Representation of a (7, 2 7 ,4) Code 



3 Classification 

Since we have shown in Theorem [T] that, except for some special cases, additive 
codes over GF(4) are equivalent to directed graph codes, it follows that to classify 
such codes, we only need to consider directed graphs. Furthermore, we only need 
to consider non-isomorphic graphs, since if two directed graphs are isomorphic, 
the corresponding directed graph codes are equivalent by a permutation of co- 
ordinates. All non-isomorphic connected directed graphs on up to 7 vertices can 
be generated in a few hours by using tools provided with the software package 
nauty [TH]. For an enumeration of these graphs, see sequence A003085 in The 
On-Line Encyclopedia of Integer Sequences [20] . 

Connected graphs correspond to indecomposable codes. A code is decompos- 
able if it can be written as the direct sum of two smaller codes. For example, let C 
be an (n, 2 n , d) code and C an (n' , 2 n ',d') code. The direct sum, C®C = {u\\v 
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u S C,v S C'}, where || means concatenation, is an (n + n',2 n+n ,mm{d,d'}) 
code. It follows that all decomposable codes of length n can be classified easily 
once all indecomposable codes of length less than n are known. The total number 
of codes of length n, i„, is easily derived from the numbers i n of indecomposable 
codes by using the Euler transform [2"T] . 

h 



To check whether two additive codes over GF(4) are equivalent, we use a 
modified version of an algorithm originally devised by Ostergard [17] for checking 
equivalence of linear codes. We show that this method also works for additive 
codes. An additive code over GF(4) is mapped to an undirected colored code 
graph in the following way. (Note that this representation is not related to the 
directed graph representation defined previously.) First, we find a set of vectors 
of some weights that generate the code. Often, the set of all vectors of minimum 
weight d will suffice, otherwise, we add all vectors of weight d + 1 , and then all 
vectors of weight d + 2, . . ., as necessary. For each vector Ci in the resulting set, 
add a vertex Vi to the code graph. Also add n sets of three vertices, where n is 
the length of the code. Each set of three vertices represent the non-zero elements 
{1, w, u! 2 } e GF(4) in one coordinate. In every set, each of the three vertices is 
connected to each of the two other by undirected edges, to form a cycle. (This 
corresponds to the fact that any permutation of the symbols {1, u>,lj 2 } in each 
coordinate of the code gives an equivalent code.) Let the vertices Vi have one 
color, and the in other vertices have a different color. Add edges between vertex 
Vi and the n 3-cycles corresponding to the codeword Cj. E.g., if Cj has u> in 
coordinate j, then there is an edge between Vi and the element labelled lu in 
the jth 3-cycle. As an example, Fig. [5] shows the case where c\ — (ui, ui, . . . , w). 
The resulting code graph is then canonized, i.e., relabelled, but with coloring 
preserved, using the nauty software [19] . If two graphs are isomorphic, their 
canonical representations are guaranteed to be the same. Hence, if two codes 
are equivalent, their canonical code graphs will be identical. Furthermore, as 
an alternative to the method described in the introduction, we can find the 
automorphism group of a code as the automorphism group of its code graph, 
i.e., the set of all vertex permutations that map the code graph to itself. 

To classify codes of length up to 7, we take all non-isomorphic connected 
directed graphs, map them to codes, and canonize the corresponding code graphs 
as described above. All duplicates are removed to obtain one representative from 
each equivalence class. The special form of the generator matrix of a directed 
graph code makes it easier to find all codewords of small weight. If C is generated 
by C — r + luI, then any codeword formed by adding i rows of C must have 
weight at least i. This means that we can find all codewords of weight i by only 



= J^did 

d I n 
= Cl 

1 / ^ 

= - C„ + Cfct„- 



k=l 
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considering sums of at most i rows of C. This property also helps when we want 
to find the minimum distance of a code. Furthermore, if we wanted to exclude 
codes with minimum distance one from our classification, it would suffice to 
exclude graphs where some vertex has outdegree zero, since this would imply 
that there is a row in the generator matrix with weight one. 

We also note another special property of directed graph codes: Given a di- 
rected graph code, C, with generator matrix r + u>I, it can be verified that the 
additive code over Z4 given by 2r + I always has the same weight distribution 
as C. We may therefore replace the elements from GF(4) with elements from Z4 
by the mapping 1— > 0, 1 t— > 2, u> t— > 1 , to 2 1— > 3. 

TableQ]gives the number of half-rate additive codes over GF(4) by length and 
minimum distance. Note that only indecomposable codes are counted, and that 
the special cases in Theorem[TJ that do not have a directed graph representation, 
are not included. A database containing one representative from each equivalence 
class is available at http://www.ii.uib.no/~larsed/directed/. Tableland 
Table [3] give the numbers of formally self-dual and isodual codes. (Note that 
the 240 formally self-dual (6, 2 6 , 3)-codes and 3 formally self-dual (7, 2 7 , 4)-codes 
were also found by Han and Kim g].) For completeness, we include in Tabled] 
the number of self-dual codes, although we have previously classified these up to 
length 12 [H]. 



4 Circulant Directed Graph Codes 



Since it is infeasible to study all half-rate additive codes of lengths much higher 
than those classified in the previous section, we restrict our search space to codes 
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Table 1: Number of Half- Rate Additive Codes over GF(4) 



d\n 2 


3 


4 


5 


6 


7 


1 1 


4 


27 


322 


8509 


686,531 


2 1 


3 


21 


262 


9653 


1,279,641 


3 




1 


9 


644 


253,635 


4 








1 


3 


Total 2 


7 


49 


593 


18,807 


2,219,810 



Table 2: Number of Formally Self-Dual Additive Codes over GF(4) 



d\n 2 


3 


4 


5 


6 


7 


1 1 


1 


10 


55 


1082 


36,129 


2 1 


2 


12 


79 


2348 


192,201 


3 




1 


5 


240 


55,711 


4 








1 


3 


Total 2 


3 


23 


139 


3671 


284,044 



Table 3: Number of Isodual Additive Codes over GF(4) 



d\n 


2 


3 


4 


5 


6 


7 


1 


1 


1 


8 


27 


344 


3243 


2 


1 


2 


10 


45 


598 


8517 


3 






1 


5 


124 


3299 


4 










1 


3 


Total 


2 


3 


19 


77 


1067 


15,062 



Table 4: Number of Self-Dual Additive Codes over GF(4) 



d\n 2 3 4 5 6 7 

2 1 1 2 3 9 22 

3 114 

4 1 



Total 1 1 2 4 11 26 
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corresponding to graphs with circulant adjacency matrices. A matrix is circulant 
if the ith row is equal to the first row, cyclically shifted i — 1 times to the right. 
The generator matrix of a directed graph code is obtained by setting all diagonal 
elements of the circulant adjacency matrix to lo. There are 2™ _1 such codes of 
length n, some of which may be equivalent. We also consider bordered circulant 
adjacency matrices: Given a length n circulant graph code with generator matrix 
C, we obtain a code of length n + 1 with generator matrix 

1 

: C 
V J 

There are 2 n_2 such codes of length n, some of which may be equivalent. For each 
n up to 26, we have counted, up to equivalence, all circulant and bordered circu- 
lant directed graph codes. The number of codes of the highest found minimum 
distance for each n is given in Table [5l A database of these codes is available at 
http : //www. ii .uib.no/~larsed/directed/ 

Proposition 2. A circulant or bordered circulant directed graph code will always 
be isodual. 

Proof. A circulant directed graph code of length n has generator matrix 



( 10 


a i 


a 2 ■ 


■ a n 


a„ 


LO 


ax ■ 


■ 0,-n-l 


O-n-l 


a„ 


LO ■ 


■ O-n-2 


\ Oi 


«2 


a 3 ■ 
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where (0,1,0,2, . . . ,a n ) is any binary sequence of length n — 1. It follows from 
Prop. [T] that the dual code is generated by 



/lo 




■ ai\ 


a\ 


lo a n ■ 


■ a 2 




a± lo 


■ a 3 




O-n-l a n-2 ■ 


■ lo) 



We can obtain C T from C by reversing the order of the columns, and then re- 
versing the order of the rows. Permuting rows has no effect on the code, and 
permuting columns produces an equivalent code. Hence the code must be equiv- 
alent to its dual. The same argument holds for bordered circulant codes, except 
that the first row and column remain fixed. □ 

If C is a circulant directed graph code of length n, then |Aut(C)| must be 
divisible by n, since the structure of the generator matrix ensures that cyclically 
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Table 5: Number of Circulant and Bordered Circulant Directed Graph Codes of High- 
est Found Minimum Distance 



n 


Max d 


# Codes 


# Self-dual 


2 


2 


1 


1 


3 


2 


2 


1 


4 


3 


1 





5 


3 


3 


1 


6 


4 


1 


1 


7 


4 


2 





8 


4 


11 


1 


9 


4 


22 


2 


10 


5 


4 





11 


5 


21 





12 


6 


2 


1 


13 


6 


2 





14 


6 


54 


3 


15 


6 


325 


3 


16 


7 


1 





17 


7 


9 


1 


18 


8 


1 


1 


19 


7 


1366 


4 


20 


8 


4 


3 


21 


8 


42 





22 


8 


1328 


17 


23 


8 


8027 


2 


24 


9 


1 





25 


9 


25 





26 


9 


1877 
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shifting all codewords will preserve the code. Similarly, if C is a bordered circulant 
directed graph code of length n, then |Aut(C)| must be divisible by n — 1, since 
the code is preserved by fixing the first coordinate and cyclically shifting the last 
n — 1 coordinates. 

With our method, we are able to find new codes, since the existence of 
near-extremal formally self-dual codes of lengths 11 and 13 was previously an 
open problem [3]. We also answer the question of the existence of an additive 
(13, 2 13 , 6) code in the positive. Parameters for optimal additive codes over GF(4) 
were determined for n < 12 by Blockhuis and Brouwer [22], and for n < 13 by 
Bierbrauer et al. [18] , by using geometric descriptions of codes. Bierbrauer et al. 
found an additive (11, 2 11 , 5) code, but posed as an open question the existence 
of an additive (13, 2 13 , 6) code. 

There are at least 21 formally self-dual (11, 2 11 , 5) codes, with generators 



available at http://www.ii.uib.no/~larsed/directed/ We find codes with 



five different weight enumerators: 

W 11A (l,y) = l + 55y 5 + 242y 6 + 275y 7 + 495y 8 + 605y 9 + 286y 10 + 89y n , 
Wii, 2 (l, y) = 1 + 66y 5 + 198y 6 + 330y 7 + 495y 8 + 550y 9 + 330y 10 + 78y n , 
Wxx, 3 (l, y) = 1 + 70y 5 + 182y 6 + 350y 7 + 495y 8 + 530y 9 + 346y 10 + 74y n , 
W X1A {1, y) = 1 + 75y 5 + 162y 6 + 375y 7 + 495y 8 + 505y 9 + 366y 10 + 69y n , 
Wu )5 (l, y) = 1 + 77y 5 + 154y 6 + 385y 7 + 495y 8 + 495y 9 + 374y 10 + 67y n . 

There are five codes with weight enumerator VFn^. Of these, three have au- 
tomorphism groups of order 11, one has 10 automorphisms, and one has 110 
automorphisms. There are four codes with Wn t 2, all with 11 automorphisms. 
There are two codes with Wua, both with 10 automorphisms. There are five 
codes with Wii t i, all with 10 automorphisms. There are five codes with ffu^, 
all with 11 automorphisms. 

We have found two formally self-dual (13, 2 13 , 6) codes. £13,1 is generated by 
all cyclic shifts of (ojIOIOOIIIOOOO) and has an automorphism group of order 13. 
^13,2 is generated by all cyclic shifts of (wlllOlllllOlO) and has an automor- 
phism group of order 78. Both these codes have the same weight enumerator: 

W 13 (l, y) = l + 247y 6 + 481y 7 + 936y 8 + 1625y 9 + 2197y 10 + 1755y n + 
715y 12 + 235y 13 . 

Note that for several lengths, there are no self-dual codes among the circulant 
and bordered circulant codes with highest minimum distance. The best known 
self-dual codes of length 24, 25, and 26 have minimum distance 8. We find a single 
isodual (24, 2 24 , 9) code generated by the cyclic shifts of (wOllOllllllllllOlOOOOHO) 
with automorphism group of order 72 and weight enumerator 

W 2i (l, y) = 1 + 1752y 9 + 8748y 10 + 26064y n + 81408y 12 + 232776y 13 + 

573516y 14 + 1119264y 15 + 1869777y 16 + 2676456y 17 + 3096804y 18 + 
2959056y 19 + 2204568y 20 + 12 5 54 1 6y 21 + 520740y 22 + 134 2 08y 23 + 
16662y 24 . 
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We also find 25 isodual (25, 2 25 , 9) codes with 25 different weight enumerators, 
and 1877 isodual (26, 2 26 , 9) codes with 1865 different weight enumerators. 

We have previously studied circulant undirected graph codes [11] . There are 
only 2 T a - 1 such codes of length n, due to the fact that the generator matrix 
must be symmetric. Gulliver and Kim [23] also performed a computer search of 
circulant self-dual additive codes over GF(4), but their search was not restricted 
to graph codes. 

A particularly interesting type of circulant code is a type of quadratic residue 
code [5]. The length of such a code must be a prime p. When p = 1 (mod 4), 
the quadratic residue code will be self-dual, and the corresponding undirected 
graph is known as a Paley graph. When p = 3 (mod 4), the code will only be 
isodual. The first row of the generator matrix of the code is (ojJi, . . . 
where ?j = 1 if i is a quadratic residue modulo p, i.e., if x 2 = i (mod p) has a 
solution x G GF(p). Otherwise, l - L = 0. Many codes with high minimum distance 
can be obtained from this construction. For instance, by bordering quadratic 
residue codes, as described above, we obtain self-dual (6, 2 6 ,4), (14, 2 14 , 6), and 
(30,2 30 ,12) codes, and isodual (4,2 4 ,3), (8,2 s , 4), and (12, 2 12 , 6) codes. 

Acknowledgement The authors would like to thank Jtirgen Bierbrauer for helpful 
comments. This research was supported by the Research Council of Norway. 
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